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Abstract 


A class of problems in air traffic management asks for a scheduling al- 
gorithm that supplies the air traffic services authority not only with a 
schedule of arrivals and departures, but also with speed advisories. Since 
advisories must be finite, a scheduling algorithm must ultimately produce 
a finite data set, hence must either start with a purely discrete model 
or involve a discretization of a continuous one. The former choice, often 
preferred for intuitive clarity, naturally leads to mixed-integer programs, 
hindering proofs of correctness and computational cost bounds (crucial 
for real-time operations). In this paper, a hybrid control system is used 
to model air traffic scheduling, capturing both the discrete and continu- 
ous aspects. This framework is applied to a class of problems, called the 
Fully Routed Nominal Problem. We prove a number of geometric results 
on feasible schedules and use these results to formulate an algorithm that 
attempts to compute a collective speed advisory, effectively finite, and 
has computational cost polynomial in the number of aircraft. This work 
is a first step toward optimization and models refined with more realistic 
detail. 
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1 Introduction 


Management of air traffic that is confined to a network of nominal routes 
requires that each aircraft be scheduled to at least a few crucial way- 
points and arrive at each waypoint as near as possible to the scheduled 
time. This requirement gives rise to an entire class of scheduling prob- 
lems, which vary by airspace configuration and various flight restrictions, 
such as speed- or altitude ranges, aircraft type, and weather. Since the 
airspace is often modeled as a directed graph, and a schedule as a chart 
that estimates for each aircraft a time of arrival at each of the finitely 
many waypoints, models of centralized scheduling of arrivals are fre- 
quently also sought in this, discrete, setting [?,?,?,?,?,?,?]. 

In such a model, aircraft are scheduled only at specific points (the 
vertices of the graph), with no time stamps maintained or imputed of 
the aircraft’s position or speed away from such a point. One weakness 
of such inherently discrete models is their principal inability to assure, 
without adding artificial constraints, pairwise separation between the 
aircraft continuously in time, as is operationally required [?]. Another 
weakness, specific to the models that lead to a mixed-integer program 
(MIP) [?], is the difficulty of finding an algorithm that performs suffi- 
ciently better than the TVP-hard worst case of a general mixed-integer 
linear program (MILP) [?], although a number of the algorithms devel- 
oped for solving a MILP modeling air traffic have been shown to perform 
with a running time polynomial in the number of aircraft (see, for ex- 
ample, [?,?]). However, while sufficient constraints can be imposed in a 
discrete model to assure separation at all times, the authors are unaware 
of published discrete models that assure such separation and show com- 
putational performance suitable for aircraft traffic in a terminal airspace. 

In this paper, we approach the scheduling of air traffic as a time- 
continuous problem. The hybrid control system (HCS) framework pro- 
posed in [?] is used for modeling air traffic motion continuously in time 
in a broad class of airspace regions represented as directed multigraphs. 
A schedule sought in this framework is, thus, a speed advisory, i.e. a 
continuous time parameterization of each aircraft’s motion along a pre- 
determined path. The pairwise separation requirements are imposed in 
continuous time. 

This framework has two main merits. First, it is capable of captur- 
ing realistically the following aspects of the research field of Air Traffic 
Management (ATM): separation requirements for every pair of aircraft at 
every instant in time, speed restrictions, off-nominal routes, and uncer- 
tainty in the actual times of an aircraft’s arrival at a waypoint. Second, 
it rests on the apparatus of HCS, which is relatively well-understood the- 
oretically [?,?,?,?,?,?,?] and furnished with general numerical methods 
(see, e.g., [?]). While no optimization is pursued in this paper, and only 
feasible trajectories are sought for a given HCS, the framework can serve 
as a setting for problems in optimal control [?]. 
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The contribution of this paper is an algorithm for the class of Fully 
Routed Nominal Problems (section ??) in ATM, where the goal is to 
guide a given finite set of aircraft, each along its pre-determined path, 
out of the route network in a separation-compliant way. The algorithm, 
formulated in section ??, constructs, or reports failure to construct, a 
feasible (i.e. , compliant with the separation and speed range constraints) 
speed advisory that guides all the aircraft as required. The computa- 
tional cost is shown to be bounded by a low-degree polynomial in the 
number of aircraft. 


2 Background 

The inherently discrete models mentioned in section ?? have historic ori- 
gins. Planning of air transportation, from the time it came into mass 
use and to the time of this writing, has in practice been based on two 
doctrines, i), that scheduling and separation assurance are two separate 
tasks and, ii), that a schedule need not synthesize for each aircraft a 4-D 
aircraft trajectory (i.e., a time-parameterized curve in 3-D space), but 
only specify an appropriate set of waypoints and the times at which they 
should be reached by the aircraft. (Current air transportation operations 
rely on Air Traffic Control (ATC) [?] to devise and issue clearance in- 
structions to each aircraft and to modify the schedule when necessary in 
order to maintain aircraft separation.) While automated trajectory syn- 
thesis has received some attention from researchers (see [?] and references 
within) and is finding its way into air traffic operations (see, e.g., [?,?]), 
doctrine i) remains a dominant practice: currently such synthesis does 
not meet, simultaneously, a given schedule and the time-continuous pair- 
wise separation requirement. 

As these conservative tendencies have persisted, increases in traf- 
fic demand have hampered ATC’s mission-to promote “safe, orderly, 
and expeditious flow of traffic” [?]-in all aspects except safety. To ad- 
dress these newly arising difficulties of ATC, researchers have sought 
to minimize various performance indices. Examples of such indices are 
makespan [?,?], fuel burn [?,?], average or total delay [?,?], excess 
separation [?], and deviation from preferred time of arrivals and depar- 
tures [?]. With this broadened focus, researchers have approached var- 
ious types of air traffic scheduling problems in terminal airspace using 
such frameworks as dynamic network flows, constraint-based scheduling, 
and scheduling combined with separation considerations; a detailed sur- 
vey, references, and discussion of this past research can be found in [?]. 

Further work on, and wider use of, automated generation of 4-D tra- 
jectories, as envisioned in the future [?], will raise the question of how 
the responsibility for scheduling and for separation should be divided 
between the automation system and the human controller. Having 4-D 
trajectory generation handled by the automation, and separation assur- 
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ance by the controller, would leave the controller with the full weight 
of separation responsibility while considerably restricting their freedom 
in the choice of aircraft maneuvers. Such a practice is likely to limit 
the benefits associated with the use of 4-D trajectories and likely to 
be unacceptable to air traffic controllers. This consideration has led us 
to consider an investigation, undertaken herein, of algorithms for 4-D 
trajectory generation that aim, as a controller does today, to maintain 
separation along the entire continuous trajectory for each aircraft. 


3 An HCS Model of an Airspace 



Figure 1. A portion of the LAX terminal airspace for arriving traffic. 
The thin arrows indicate the traffic directions. 

For the operational purposes of ATM, an airspace is regarded as a region 
with a set of marked points, called waypoints. Some pairs of waypoints 
are connected by route segments. (For example, a portion of the airspace 
surrounding the Los Angeles International Airport is shown in Fig. ??.) 

The same pair of points may be connected by more than one segment. 
In practice, a route segment directly connecting a pair of waypoints is 
rectilinear, but this assumption is made here only in section ??. This 
setting is naturally modeled by a directed multigraph [?]: the waypoints 
are the vertices, while the route segments are the edges. A sequence of 
route segments with every consecutive pair adjacent is called a route and 
is none other than the graph-theoretic path [?]. 

Although the application central to this paper is in ATM, the follow- 
ing construct is, in principle, suitable for any finite set of moving agents 
whose movement is confined to a route network. Assume one is given 


5 


the following data: 


1. A directed multigraph G — (V, E ), each vertex v G V being a point 
in a Euclidean space E of dimension 2 or 3. If e G E is an edge from 
v\ G V to V 2 G V, then the nominal route segment from waypoint 
v\ to waypoint V 2 is a curve in E, connecting v\ to All such 
curves will henceforth be assumed rectifiable [?, section 4.6-9] and 
capable of a parameterization which is continuous and piecewise 
continuously differentiable. A cusp in the curve can be traversed 
with the assumption (made throughout this paper) that inertia is 
neglected, and approximately smoothed if inertia is included. The 
out degree and indegree of a vertex are allowed arbitrary values. To 
emphasize the aerospace context, the multigraph G = (V,E) will 
be called a route network. 

2. A finite set A — {1, . . . , A}. of moving agents a G A in G. If agent 
a is moving along edge e, the agent’s position is specified by the 
arc length coordinate along a path containing e; the choice of 
arc length coordinates will be clarified as the model is developed 
and will always have the coordinate increase along the direction of 
the edge. 

3. For each agent a G A, a specification of the agent’s initial position, 
which is a point in G specified, for example, by an edge in G and 
a fractional distance along that edge. 

4. The inertia- free dynamical law [?] (henceforth the dot denotes dif- 
ferentiation with respect to physical time t) = s^ , where the 

’ s are the corresponding speeds, describing the motion of those 
agents a that have not yet exited the system. In what follows, and 
with the details provided in the following paragraph, the coordi- 
nates Xg will play the role of state variables ; the speeds , of the 
control variables. 

5. State constraints: the separation requirement for each pair of agents. 
This requirement is described mathematically, in terms of the co- 
ordinates x a , in section ??. A state in which at least one pair of 
agents violates the separation requirement will be called a conflict- 
ing state. 

6. Control constraints: bounds on the speeds . 
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(A) (B) (C) 



Figure 2. Examples of three discrete modes with edge set E = 
{ei, e 2 , . . . , e^} and moving agent set A = {1, 2, 3} ( A = 3). The agents 
are shown as numbered gray squares. The discrete modes shown are: (A) 
Mi : Mi(l) = ei, Mi( 2 ) = ei, Mi (3) = e 4 ; (B) M 2 : M 2 (l) = e 4 , M 2 ( 2 ) = 
e 2 , M 2 (3) = e 4 ; (C) M3 : M3(l) = ei, M3( 2 ) = e 3 , M3 (3) = e 4 . 

One choice of a hybrid control system (HCS) [?] fit for modeling 
coordinated multi-agent motion in a route network G was introduced 
in [?]. In this HCS, one defines each discrete mode (in [?], location) as 
a mapping fi : A — » E, which specifies for each agent a the edge /i(a) 
currently occupied by that agent. Examples of three discrete modes on 
a 6-edge route network with 3 moving agents are shown in Fig. ??. 

(A) (B) (C) 



Figure 3. (A-C) The continuous state spaces correspond- 

ing to discrete modes /ii, /i 2 , /i 3 from Fig. ??. Among the discrete modes 
the system may enter from (i\ are fi 2 and /i 3 , accordingly as agent 2 en- 
ters e 2 or e 3 . Conflicting states not shown. (D) A gluing of to X^ 2 
and to X /13 . Glued, all three state spaces share a face, but no two have 
any other points in common. Conflicting states not shown. 

Denoting the arc length parameterization of edge /a(a) by the 
dynamical law corresponding to discrete mode /i is given by the system 
x* = a G A , of differential control equations. Suppose each x “ 

ranges from 0 to the length l(/a(a)) of edge /i(a). The set X M of all 
continuous states = (x^) ae ^ = x 2 ^ . . . , x^j in discrete mode 

fi is then obtained by taking the A-dimensional parallelotope given by 
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the Cartesian product Y[ a [0,l(g(a))] an d removing from it all states 
that violate the separation constraint for at least one pair of agents 
(conflicting states). (The set A f is an example of a roadmap coordination 
space , defined in [?, section 1.2].) For example, Fig. ??A shows the A Ml 
for discrete mode gi from Fig. ??A; the conflicting states are not shown 
in Fig. ??. 


One must then specify, for each pair of states G A^, xy G Ay, 
whether it is possible for the system to go from one to the other, thus 
switching from one discrete mode to another, instantly. The set, denoted 
S^y, of all such state pairs (x M ,xy) for a given pair of discrete modes 
g,g! is called a switching set [?]. For this HCS, the switching sets are 
completely determined by the following condition: (x M ,xy) G S^y if 
and only if there is at least one agent, a, that in discrete mode pi occupies 
a vertex v G V with coordinate x “ = l ( 11 (a)), and in discrete mode g! 
occupies that same vertex with coordinate x^, = 0. Thus, as the system 
switches from g to g', agent a exits edge g(a) and enters edge g'(a). 
Geometrically, the switching sets are described by a topological gluing 
[?] (similar to that used in [?, section 3]) between the corresponding 
continuous state spaces X M and Ay: the right endpoint of [0, l(g(a))] 
is identified with the left endpoint of [0, l(g'(a))]. If the transition from 
g to g! involves a change of edge for more than one agent, then each of 
these agents effects a gluing. Thus, the face x^-^ = l(g(a)) of A M is 
to be glued to the face = 0 of Ay for each discrete mode g! that 

allows a direct transition from g. (For example, in Fig. ??(A-C), the 
“top” (in the sense of the Figure) face of A m is glued to the bottom face 
of A^ 2 and of A M3 ; see Fig. ??D.) The state space obtained by carrying 
out all such gluings, and denoted A, is called the global continuous state 

glue 

space. Denoting the gluing operation by the non-standard symbol U , 
glue 

one has X — U M A^. Since the set of the conflicting states in each 
A m is closed [?], as it will be defined below using sharp inequalities, 
two states x M ,xy identified as a result of the above gluing are either 
both conflicting or both separation-compliant. It follows from the above 
that every global continuous state space is a chain of dimension A on 
R/ 4 [?, section 35E] with the conflicting states removed, but generally 
fails to be a differentiable manifold [?]. 


The only discrete control variable appearing in this HCS would thus 
be the one indicating, for an agent reaching a vertex of outdegree > 2, 
which edge the agent is to take (or whether the agent should exit the 
system) when leaving that vertex. 
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The Geometry of a Minimal Separation Re- 
quirement 



Figure 4 . Agents 1,2 on their respective rectilinear edges ei,e2, which 
share a common vertex, taken as the origin 0 in R 2 . The orientation 
of the edges is not specified. (A) The unit vectors ai,a2 are collinear 
with the respective edges, but their directions do not necessarily agree 
with the edges’ orientations. (B) With suitably chosen scalar coefficients 
ci, C2, the vectors ciai and C2a2 are the respective position vectors of the 
two agents. 

In some types of transportation, including aircraft and trains, every pair 
of moving agents must be separated by a distance no smaller than a pre- 
determined minimal separation. Such requirements can be asymmetric 
and anisotropic; for a detailed discussion, see [?, section II]. Here we 
will characterize states that violate the separation requirement using 
conservative approximations, addressing only the following asymmetry: 
if two moving agents are in-trail (i.e., one is directly following the other 
along a route segment which is not necessarily in a horizontal plane) and 
have different weight classes, then the minimal separation can depend on 
the weight class of the agent in front. (For example, if an aircraft of type 
Heavy is directly following one of type Small, the required separation is 
smaller than if the aircraft were interchanged.) To capture this potential 
asymmetry, for each pair aq,<a2 of agents with the first one leading, we 
introduce the minimal separation r ai ^ a2 . If the asymmetry takes place, 
it can be written 

T‘a±,a2 7*“ ^0-2 ,Qq (i) 

We now calculate the set of all conflicting states in a discrete mode /i of an 
HCS described above; recall that these are the states where at least two 
agents violate the separation requirement. The scenario shown in Fig. 
??A has two agents on two different rectilinear edges, which need not he 
in a horizontal plane, with a common vertex and no specified orientation. 
(If the edges are curvilinear with low curvature near a common vertex 
or intersection, these portions can be approximated by linear segments; 
otherwise, the analysis becomes considerably more complicated.) 
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Remark 4.1. Since edge orientation is not specified, Fig. ?? describes 
four cases : both agents are moving toward the common vertex, both mov- 
ing away from the common vertex, and two more cases in which one agent 
moves toward, and the other away from, the common vertex. 

(A) (B) 



Figure 5. (A) An example of two elliptical sectors in the ciC 2 -plane 

corresponding to conflicting states. (B) An example of two stripes in the 
ciC 2 -plane corresponding to conflicting states of two agents on the same 
edge. 

We will use the Euclidean inner product [?] (•, •) and the correspond- 
ing norm [?] || • || in the 2-D space containing the two edges. Pick the 
common vertex as the origin and the unit vectors ai , a 2 as the basis 
vectors that, regardless of the edge orientations , point from the origin 
toward the respective agents. With suitable scalars ci,C 2 , the vectors 
ciai and C 2 a 2 are the agents’ respective position vectors. The squared 
distance between the two agents, denoted by D 2 , is 

D 2 = ||ciai - c 2 a 2 || 2 = (ci) 2 + (c 2 ) 2 - 2cic 2 (ai,a 2 ) 

Equating D 2 to the squared minimal separation, say, r\ 2 , we obtain the 
equation 



of an ellipse in the ciC 2 -plane. The corresponding set of conflicting sets 
is described by the elliptical sector obtained by intersecting the open 
octant C 2 > ci > 0 with the ellipse-bound region D 2 < r\ 2 (respectively, 
D < r\ fi) when agent 1 (respectively, agent 2) is the one closer to the 
origin. The role of the angle 6 between the edges ei, e 2 in both sectors is 
the equality (ai,a 2 ) = cos(0). An example of two such sectors is shown 
in Fig. ??A. The asymmetry of the gray-shaded region about the dashed 
diagonal is the asymmetry (??). 

In each of the four cases listed in Remark ??, the respective contin- 
uous state coordinates x^x^ of agents 1,2 in discrete mode fi map to 
the coefficients ci,C 2 , as follows: 

1. If both agents are moving toward the common vertex, then x ™ = 
l(e a ) — c a. for a — 1, 2. 

2. If both agents are moving away from the common vertex, then 
xfif — c a for a — 1, 2. 
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3. If agent 1 is approaching, and agent 2 going away from, the common 
vertex, then x * = Z(e i) — ci,x^ = C 2 - 

4. If agent 2 is approaching, and agent 1 going away from, the common 
vertex, then x 1 — ci,a^ = l(e 2 ) — C 2 . 

If 6 > 90°, then in the last two cases fi allows only one in-trail sequence, 
so the minimal separations used for the two sectors in Fig. ??A are 
equal. If the two agents 1,2 are on the same edge, then the set in 
the ciC 2 -plane of the conflicting states appears as in Fig. ??B and is, 
in accordance with (??), asymmetric about the dashed diagonal. The 
mapping from the continuous state coordinates x^x 2 ^ to the coefficients 
ci,C 2 is constructed analogously to the above four cases. 


(A) 



(B) 





Figure 6. An example of two agents whose paths overlap. The black 
star shows the beginning of the overlap in (A) and the corresponding 
state (both agents being at that point) in (B); the white star, the end 
of the overlap in (A) and the corresponding state (both agents being at 
that point) in (B). The system, shown in (A), has seven discrete modes 
with both agents in the transportation network. Each mode’s set of 
separation-violating states, shown in (B) as a connected [?] gray region, 
is “glued” to some of the others. The result of the gluing is the connected 
region shown in (C). 


The above calculation is illustrated, for an example of two moving 
agents, in Fig. ??. Each discrete mode’s set of conflicting states is 
shown as a connected [?] gray region. For dimension A above 2, one must 
compute for each pair of agents the set of states violating the separation 
requirements. Each such set is a cylinder, or union of cylinders, with the 
base shaped as shown in Fig. ??C, in the total state space U^X^. We 
note that the set of all separation- violating states in is cylindrical 

in the sense of [?, Definition 2.2], the latter definition a key requirement 
for the applicability of a number of theoretical results of [?]. 


5 An Alternative: Path-Based Transportation 
HCS 

The HCS defined in this section will be instrumental in a precise formu- 
lation of the Fully Routed Nominal Arrival problem, defined in section 
??. Assume the data ??)-??), listed in section ??, are given. 
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1. For each agent a G A, let e INIT ’ a denote the edge occupied initially 
by agent a. (Two or more agents can occupy the same edge.) 

2. Let V (e INIT '^ be the se t Q f a p paths in the multigraph that 

begin with the edge e INIT]a and end with an edge that contains 
an acceptable destination for a. Denote the length of a path p G 
V ^e INIT i a ^J by l(p). 

3. Define each discrete mode as a mapping that assigns each agent 

to a path in V ^ e INIT ^ m l n more detail, v is a mapping from 

the set A of moving agents to the union U a V ^ e INIT '^ such that 
v(a) G V (e INIT '°^j for each a G A. 

4. In each z/, have the arc length coordinate evolve according to 
the dynamical law 

Vv = ote A, ( 2 ) 

where s “ is the control variable corresponding to the agent’s speed 
of motion along the path. 


5. For each agent a and each discrete mode /i, specify the speed con- 
straint 


ga ; mm ^ $ a ^ ^cc;max 


(3) 


where 5"^ ;max ] is the speed range allowed for agent a in 

discrete mode fi. 


Definition 5.1. The HCS defined in this section will be called a path- 
based transportation HCS. 

In a path-based transportation HCS, every discrete mode v corre- 
sponds to a possible routing of the agents in A , and every execution 
occurs in only one discrete mode. It follows that the discrete controls 
in a path-based transportation HCS are, effectively, identified with the 
discrete modes, and play but a trivial role. 


6 A Geometric Theory of the Fully Routed Nom- 
inal Arrival Problem 

The geometric results obtained in this section will be used in the formu- 
lation of the algorithm in section ?? to prove correctness and bounds on 
computational costs for some of the steps. 


12 



6.1 Problem Definition and Geometric Properties 

In this section, we describe and analyze theoretically a subclass of the 
HCS defined above, the analysis being a foundation for the algorithm in 
section ??. This subclass represents a most commonly modeled scenario 
of air traffic, which is arriving in merging and, possibly, diverging flows, 
with no multiple routing options: each arriving aircraft follows without 
deviation a prescribed path in the route network. Consequently, the 
resulting path-based transportation HCS (definition ??) has only one 
discrete mode. Henceforth, the moving agents are aircraft. The subscript 
v will be dropped for brevity. 

Remark 6.1. The assumption that each aircraft is assigned a route may 
not always hold. In cases when such assignment (known in Air Traffic 
Operations as routing,) is not given, but is sought as part of the problem, 
the corresponding HCS has multiple discrete modes. The theory and 
algorithm developed in this section can then be used for each of those 
discrete modes, and the result computed in one of the modes selected 
later as final. See [?] for an analysis of computational costs and the 
possibility of parallel computation of routing. 

The problem consists in computing a control strategy (s a (t)) ae ^ 4 that 
will, in accordance with dynamical law (??), take each aircraft to the end 
of its prescribed path. To ATM operations, such a control strategy would 
serve as a speed advisory. This problem will be called the Fully Routed 
Nominal Arrival problem. It consists in finding feasible trajectories. No 
optimization is pursued. 

Section ?? gives a list of assumptions made to simplify the compu- 
tations. As is shown below, these assumptions allow a solution by an 
algorithm that is noniterative and proven to halt in polynomial time (in 
the number of aircraft). Because of the urgency that can accompany 
ATM decisions, this bound on computation time is a highly desired fea- 
ture of every solution candidate to be deployed as an operational tool in 
the field, especially for automated tools responsible for both separation 
assurance and scheduling. 

6.2 Simplifying Assumptions and Their Implications 

Assumption 6.1. All control strategies considered before section ?? will 
be assumed executable with certainty, i.e. without execution error that 
may be caused by wind, transportation performance, or human factors. 
In particular, under this assumption, the aircraft is able to follow its 
prescribed route exactly, without deviation. 

An approach to including execution uncertainty in the model is dis- 
cussed in section ??. 
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Assumption 6.2. All aircraft have the same cruising speed range , 

j gmin ^maxj 

Assumption 6.3. All aircraft pairs have the same required minimal 
separation distance, henceforth denoted r. 

Assumption 6.4. During the execution of the algorithm described in 
section ??, below, no new aircraft enter the system. 

6.3 Geometry of the Global Continuous State Space 

6.3.1 The Global Coordinate System 

Since the system is fully routed, each aircraft a is assigned a uniquely 
determined path, whose length will henceforth be denoted M a , a — 
1,...,A Accordingly, the global state space Y for the problem can 
be obtained by removing from the parallelotope YlaeA^^^] all the 
conflicting states. We choose the coordinate system, for the whole of 
Y, so that the state variable y a indicates the arc length coordinate of 
aircraft a on its path and ranges from 0 to M a . Unless said otherwise, 
all state variables are regarded as coordinates referred to the standard 
basis ei = (1, 0, . . . , 0, 0), . . . , = (0, 0, ... , 0, 1). 

6.3.2 The Target Set and the Distal Boundary 



Figure 7. The case of A = 3 aircraft: (A) The distal edges and the target 
set, and (B) an illustration of the distal boundary. 
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The vertex M = (Mi, . . . ,M^) of Y will be called the distal vertex of 
Y. The faces of Y containing the distal vertex will be called distal faces. 
Distal faces of dimension 1 will be called distal edges. Since the distal 
edges correspond, one-to-one, to states in which all but one aircraft have 
exited the airspace, we define the target set of the Fully Routed Nominal 
Arrival problem as the union of all distal edges. Fig. ??A illustrates the 
distal edges for 3 aircraft. 

The union of all the ( A — l)-dimensional faces adjacent to the distal 
vertex will be called the distal boundary. An illustration of the distal 
boundary for the case of 3 aircraft is given in Fig. ??B. 


6.3.3 Polygonally Approximated Pairwise Conflict Zones: the 
Geometric Properties Used in the Algorithm 


Each pairwise conflict zone will be approximated here by a polygon that 
is the intersection of three half-planes. This approximation achieves 
considerable computational advantage, but likely incurs a cost if an op- 
timization were pursued: by enlarging the conflict zone, one restricts the 
class of considered candidate solutions. This loss would be most severe 
if at least two of the paths of the agents were to cross (the relevant 
operational concepts are same courses , crossing courses , and opposite 
courses ; see [?, section 1-2-2].) In the framework of the above model, 
such a crossing can be viewed as a short overlap (compare to Fig. ??). 
To analyze the geometry of the conflicting state sets that result from this 
approximation in the global state space X, it will be convenient to in- 
troduce the technical terms of protrusion , safe wedge , innermost conflict 
edge , and distal innermost conflict-free point. 

In the A-dimensional Euclidean space R/ 4 , define the unit vector 

d = •••>!)• 

Let H a i’ a2; -i and H a 1,a2;+1 denote, respectively, the half-spaces con- 
sisting of all states y = (y 1 , . . . , y A ) satisfying y ai — y a2 > — r ai ’ a2 and 
y°L\ _ yO .2 < r ai,a 2 ^ where r «i ,«2 j s the sum Q f the required minimal sep- 
aration distance between aircraft oq and 0 L 2 and a “buffer,” which is a 
value sufficiently large to make sure that the obtained polygonal approx- 
imation contains all of the exact pairwise conflict zone depicted in Fig. 

??C. 
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Figure 8. The shaded regions are the orthogonal projections of the half- 
spaces (A) fjai,a 2 ]+l ^ anc [ (Q) F ai,a2 onto the y ai y a2 ~ 

plane. Their intersection is the polygonal approximation, shown in Fig. 
??D, of the conflict zone for oq, o^- (D) An approximation of a pairwise 
conflict zone (Fig. ??C) by a polygonal region equal to the intersection 
of three half-planes. (E) The geometry of the safe wedges (shaded) and 
DICPs. 

Define the half-space F ai ,a2 to consist of all states y satisfying (y — 
M, — d) < where f ai ^ a2 is the orthogonal distance from M to the 

bounding hyperplane of F ai,a2 (Fig. ??). 

The pairwise conflict zone for the aircraft pair aq,02 is, then, ap- 
proximated by the open [?] polyhedral region 

1,«2 _ PJOL 1 ,^ 2 ; — 1 p p| pai,a2 

depicted in Fig. ??D. The distance f ai > a2 will be called the protrusion 
of the conflict zone (??). 

6.3.4 Safe Wedges, Innermost Conflict Edges, and Distal In- 
nermost Conflict-Free Points (DICPs) 

The intersection H ai ^ 2 '~ l D H ai ^ a 25+1 is a conflicting slab of the form 
\y ai — y a2 \ < r. By “removing” from the A-dimensional Euclidean space 
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H A all the conflicting slabs (i.e., the slabs corresponding to all pairs of 
aircraft), one obtains a pairwise disjoint collection of A\ polyhedral re- 
gions, called safe wedges. The intersection of a safe wedge with the global 
continuous state space is a conflict-free region and has a 1-dimensional 
face (a line) parallel to d. This line will be called an innermost conflict- 
free edge. The intersection of this edge with a distal face of dimension 
(. A — 1) is a point called a distal innermost conflict-free point (DICP). 
Safety wedges and DICPs for the case of A = 2 aircraft are illustrated 
in Fig. ??E. 

Consider “the positive A-dimensional quadrant” 

R >o = {(y\ y 2 , ■ ■ ■ , y A ) ■ y a > ° for ail 

of the A-dimensional Euclidean space RA For each 
the indices 1, . . . , A, the set of all y = (y 1 , y 2 , . . . , y A ) 

0 < y a{1) < y a{2) <...< yAA 

determines a subset Y a of the “quadrant” R> 0 . 

Definition 6.1. The subset Y a (in fact , the interior of a polyhedral 
angle) will henceforth be called the basic polyhedral angle corresponding 
to the permutation a. 

Note: Henceforth, the end of a mathematical proof is marked by an 
open square. 

Lemma 6.1. If the permutation a is such that the safe wedge W consists 
of all y = (y 1 , y 2 , . . . , y A ) with 

y CT (“) < r + y <T ( a+1 ) for a = 1, 2, . . . , A — 1, (6) 

then W is contained in Y a . 

Proof. Condition (??) follows from (??). □ 

Theorem 6.1. If all the required minimal distances for in-trail separa- 
tion have the same value, r, then the coordinate A-tuple of every DICP 
is the sum of M and a permutation of 0, — r, —2 r, . . . , — (A — 1 )r, and, 
conversely, every such permutation gives a DICP. 

Proof. These H-tuples are exactly those at which all constraints (??) are 
active. □ 

6.3.5 The Cone of Attainability, Feasible Trajectories, and 
Feasible Solutions 

The flow of the traffic consisting of A aircraft a = 1, . . . , A can be seen 
as a trajectory (y 1 (t), y 2 (t), . . . ,y A {t )) of the inertia-free dynamical law 

(??)• ' 


a} 

permutation a of 
such that 

(5) 
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Remark 6.2. Throughout this paper , the term collective trajectory (or, 
simply, trajectory^) will mean a time-parameterized curve 

in the global continuous state space X that satisfies ( 11 ) for some choice 
of a control strategy s (t) — (s l (t ), . . . , s A (t)). This usage of the term 
trajectory differs from that commonly seen in Air Traffic Management. 



Figure 9. An example with two aircraft. The conflict zone is as in Fig. 
11 D. (A) A collective state trajectory y (£), with initial state y(0) = yo. 

(B) The cone of attainability C(yo) (shown in the lighter shade) at yo- 

(C) The cone is positioned so that there are no feasible solutions. 


Recall assumption 11 . A segment yo y in X is said to be attainable 
if all aircraft pairs satisfy 

cmin omax 

- vo 2 ) < y ai - vo 1 < ^r(y“ 2 - Vo 2 ) (V 

A ray emanating from yo is said to be attainable if for every point y on 
the ray the segment yoy is attainable. The union of all attainable rays 
emanating from yo is a polyhedral cone, called the cone of attainability at 
yo and denoted C(yo). This cone is a pointed cone with vertex yo- A col- 
lective trajectory y(t) with initial state y(0) = yo is said to be attainable 
if it is piecewise differentiable and, at every point y of differentiability, 
the corresponding tangent vector lies in the cone of attainability at y. 
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A conflict-free attainable trajectory is said to be feasible. A 2-aircraft 
example is shown notionally in Fig. ??. 

Definition 6.2. A feasible trajectory whose final state is in the target 
set is called a feasible solution. 

The following lemma characterizes all the collective states from which 
it is safe to have all aircraft fly at the same speed. 

Lemma 6.2. If yo is a DICP, then there exists a feasible solution with 
initial state yo- Consequently, a feasible trajectory with the final state a 
DICP can be extended to a feasible solution. 

Let | • |oo denote the max norm, defined for a vector £ = (£1, £2, • • • > £a) 
in R a by |£|oo = max a |£ a |. 

Proof. The collective trajectory y (t) = yo + t.S' rnaX ]^j— d. with the a-th 
component excluded once aircraft a exits the global continuous state 
space, is conflict-free since yo lies in a safe wedge and since y is co- 
direct ional with d. □ 

6.3.6 The Geometry of FCFS: the Closest DICP to y Corre- 
sponds to the Safe Wedge Containing y 

The key result of this section, theorem ??, tells how to find (in time 
0 (A 2 )) a DICP that requires “minimal position shifting” to be performed 
on the current collective state. 

Lemma 6.3. If a = (01,02) and b = (61,62) are vectors in the first oc- 
tant of the 2 - dimensional Euclidean space R 2 with Cartesian coordinates 
yi, y2, i-e. satisfy 0 < a\ < <22 and 0 < 61 < 62, then 


aibi + a 2 6 2 > ai6 2 + a 2 6i (8) 

Proof. Both a and b lie in the positive quadrant of the plane and below 
the diagonal y\ — y2- Therefore, the angle between a and b is smaller 
than that between a and b' = (62,61) (the latter point being above the 
diagonal). Since all angles are between 0 and n/ 2 , one has (a, b) > 
(a, b 7 ), which is inequality (??). □ 

Lemma 6 . 4 . Suppose a = (aq, . . . , a a) and b = (61, ... , 6^) are vec- 
tors, with all standard components distinct, in the A-dimensional Eu- 
clidean space R a such that 0 < a\ < 02 < . . . < a a and 0 < 

b a for ol — 1 , , A. For each permutation a of { 1 , 2 , ... , A}, let ah = 

(6<t(i) 5 6^(2) , . . . , 6 ^^)). Then, of the products (a, crb) corresponding to all 
the different permutations a, a largest one is delivered by that permuta- 
tion cr* for which 6^(1) < 6^(2) < • • • < ^(A)* 
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Proof. Applying BubbleSort [?] to the numerical sequence 


h (9) 

one obtains a sequence o\, 02 , . . . , crj_i, aj — cr* of permutations of 
{1,2,..., A}, every two consecutive ones differing by a transposition 
that puts two elements of (??) in increasing order. By lemma ??, 
(a, aj- ib) < (a, ajb) for j = 2, 3, . . . , J. Thus, the largest of the prod- 
ucts (a, ajb) is given by aj = a*. □ 


The following helps locate a DICP closest to a given y. 

Theorem 6.2. If y and a DICP p lie in the same basic polyhedral angle , 
then p is a closest DICP to y. 

Proof. For this proof, change coordinates so that M is now the origin, 
and d points along the positive diagonal. If q were a DICP closer to 
y than p is, then the permutation a such that q = crp would satisfy 
(<rp,y) > (p,y), contradicting lemma ??. □ 


6.3.7 Existence and Nonexistence of Feasible Solutions 

Lemma 6.5. (A sufficient condition for the existence of a feasible solu- 
tion.) If yo is in a safe wedge, then there exists a feasible solution with 
initial state yo. 

Proof. Consider the trajectory y (t) — yo + tsy^|— d, with t increasing to 
the value T for which y(T) lies on the distal boundary. By construction, 
the trajectory is feasible, and y(T) lies in a safe wedge. The state yx, 
taken to be the initial state for a new, ( A — l)-dimensional, problem, lies 
in a safe wedge of the new ambient (A — l)-dimensional space. (For, in 
the state yx, at least one aircraft, say o-th, has exited the airspace, so 
all the pairwise conflict zones involving a are removed.) Consequently, 
the lemma holds by induction on A. □ 


Theorem 6.3. (A sufficient condition for the absence of feasible solu- 
tions.) If a pairwise conflict zone (??) (not necessarily with a largest 
protrusion) is such that C(yo)ndF ai,a2 is contained entirely in the con- 
flicting slab H ai ^ 2 '~ l p F[ a i,« 2 ;+i (Fig. ?? C), then there is no feasible 
solution with initial state yo- 

Proof. In a cone C(yo) so positioned, every attainable collective trajec- 
tory reaches a conflicting state before reaching the target set. □ 
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7 A Polynomial-Time Speed Control Algorithm 
for Computing a Feasible Solution or Report- 
ing That None Was Found 

Keep the notation introduced above. In the following algorithm, the 
computational cost (abbreviated C.c.) of each step is given in square 
brackets. 

Given: An initial state yo, global continuous state space (completely 
determined by (M — yo)), and the pairwise conflict zones (??). 

1. Find all conflicting slabs that contain yo. [C.c.: (^) = 0(A 2 ).] 

2. If yo lies in no conflicting slab, i.e., lies in a safe wedge, exit the 
algorithm, returning the speed advisory constructed in lemma ??. 
[C.c.: 0(A).} 

3. Test whether the collective state yo is in any of the pairwise conflict 
zones (??). [C.c.: 0(A 2 ).} If it is, go to step ??. 

4. Find a DICP p closest to yo- [By theorem ??, the computational 
cost of this step is 0(A log 2 A).] 

5. Compute the state yi = p + ad that lies on the boundary of 
C(y 0 ). [C.c.: O(A).} 

Remark 7.1. By construction, the state yi lies in a safe wedge, 
hence the segment yip is a feasible trajectory (lemma ?? ). If seg- 
ment yoyi is feasible, then the trajectory yoyiP is feasible and, by 
lemma ?? ; extends to a feasible solution. 

6. If segment yoyi is feasible, exit the algorithm, returning the speed 
advisory for the trajectory described in remark ??. [C.c.: 0(A 2 ). 
Computations are required to check the feasibility of yo yi .] 

7. (In this step, let Y denote the hyperplane orthogonal to d and 
passing through the origin. Orthogonal projection onto Y will also 
be denoted by Y.) Segment yo yi has been found infeasible. Use a 
fixed number of iterations of the Bisection Method to look for the 
parameter (3 that determines the point q = /3Y p such that the 
point yi;*, defined as the intersection between the line through q 
parallel to d and the boundary of C(yo), yields a feasible segment 
yo yi ; *. If this segment is found, return a speed advisory for the 
trajectory yoyi ; *P (analogously to step ??). Otherwise, go to step 
?? 

8. Exit the algorithm with the message, “Cannot produce a conflict- 
free speed advisory by speed control only.” 
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Table 1. The output format of an advisory segment. 


a/cl speed 


a/c 2 speed 


a/c A speed 


TTF 


8 Sample Numerical Solutions 

If a feasible speed advisory is found by the algorithm of section ?? suc- 
cessfully, it is output as one or more rectilinear segments, each called an 
advisory segment and conforming to the format shown in Table ??. Here 
TTF stands for the time to fly (with that speed) and indicates the time 
duration in which the advisory segment is to be traversed. However, in 
the numerical results presented below (Table ??), instead of showing all 
the individual speeds, the following abbreviated format will be used: 

(10) 


of the A speeds max. of the A speeds 


TTF 


# aircraft: 23; actual speed min & max: [21 1 .84, 250] 



nmi 


Figure 10. The initial state of a traffic scenario for 23 aircraft. To capture 
airspace adequately, an aspect ratio different from 1 is used. The thin 
arrows indicate the directions of the traffic. 

The above algorithm was applied to a number of arriving traffic sce- 
narios in the portion of the LAX airspace shown in Fig. ??. Although all 
the paths in this airspace portions converge to one runway, such conver- 
gence is not a requirement: the algorithm is applicable to all situations 
where the aircraft follow pre-determined paths. The solutions found for 
two of the scenarios, each with two advisory segments, are summarized 
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Table 2. Each of the Tables gives a 2-segment speed advisory and TTF 
for an arriving traffic scenarios. In columns 4-6 of each Table, each line 
fits format (??). 


# 

yo 

ad vis. 

min. 

max. 

seg. 

a/c 


seg. 

speed 

(kts) 

speed 

(kts) 

TTF 

(hr) 

23 

Fig. ?? 

i 

211.84 

250.00 

0.56 

2 

250.00 

250.00 

0.41 


# 

yo 

ad vis. 

min. 

max. 

seg. 

a/c 


seg. 

speed 

speed 

TTF 




(kts) 

(kts) 

(hr) 

24 

not 

shown 

i 

199.11 

250.00 

0.59 

2 

250.00 

250.00 

0.44 


in Table ??. The initial states are conflict-free. One of the two scenarios 
is depicted in Fig. ??, with circles of radius r centered at the aircraft. 


9 Discussion 

The above HCS framework addresses simultaneously the following issues 
in Air Traffic Management and adjacent fields: scheduling, speed advi- 
sory generation, conflicts (avoidance, detection, and recovery), compu- 
tational efficiency, and algorithm correctness. The theory and algorithm 
given above are mainly a proof of concept, intended to demonstrate how 
HCS can be used to obtain provably correct solutions efficiently. The 
following are some directions for further research. 

Control equations with inertia: A model more realistic than the dy- 
namical law (??) should include inertia, thus making the agents’s accel- 
eration (and not speed) the quantity directly controllable. 

Allowing new aircraft to enter the airspace during the computation: A 
reasonable direction is a search for geometric insight into those entrances 
for which the original, A-dimensional, problem has a feasible solution 
but the newly obtained one does not. The problem above, together 
with the Scheduled Routing Problem formulated in [?], suggest that 
such insight would be valuable not only in air traffic management and in 
airspace design, but also in the general field of multi-agent coordination. 
In published literature, entrance of new aircraft into the system has been 
considered in the context of the dynamic Aircraft Scheduling (Landing) 
Problem (see [?]), which, however, does not consider separation except 
at the runway. 

Conflict recovery: Recovering from a conflicting state y c in the above 
model is equivalent to finding-or establishing the absence of-a shortest 
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(or otherwise optimal) attainable trajectory from y c to a nearest safe 
wedge. Theorem ?? is instrumental here, as a nearest safe wedge is with 
a nearest DICP, and the latter is computable in time polynomial in A. 

Path control (in ATM operational terminology, Vectors): In opera- 
tional ATM practice, an aircraft is often allowed, and even instructed, 
to deviate from a given route. Such deviations are known in Air Traffic 
Operations as path control This deformation of the given route results 
in a deformation of the global continuous state space and of the pairwise 
conflict zones. Insight into the relation between the two deformations 
can be instrumental in extending the model to include path control. 

Modeling uncertainty of control execution: Control of aircraft is af- 
fected by such sources of uncertainty as human factors, aircraft perfor- 
mance, and weather. In the context of the above HCS, such uncertainty 
generally leads to stochastic dynamical law [?, chapter 14]. Among the 
most common sources of random perturbation to an aircraft’s execu- 
tion of control is wind. In the special case when each aircraft is cleared 
for a procedure that prescribes a specific route consisting of rectilinear 
segments, connected by fixed-radius turns (e.g., RF-leg types [?]), the 
aircraft’s navigation system will keep the aircraft on the route even in 
the presence of wind. Therefore, the wind affects only those components 
of the aircraft’s speed and acceleration tangential to the route. This 
suggests a stochastic dynamical law, for the Fully Routed Problem, that 
includes inertia (i.e., the control variable is the acceleration), governs 
the arc length coordinate y a and the arc length speed s a as stochastic 
processes serving as the state variables, and has the form 

dy a — s a dt , ds a = a a dt + a w dt + dW a , a G A 

where a w is the (possibly time-dependent) mean of the tangential com- 
ponent of the acceleration due to wind, and IF is a Wiener process that 
characterizes the fluctuations of the acceleration due to wind about the 
mean. 

10 Acknowledgment 

D. Davis’s research was supported by NSF Grant DGE-0707424. The 
authors thank D. Denery (NASA ARC) for suggesting the use of hybrid 
systems for ATM, W. Haskell (USC) for editing early drafts of the paper, 
H. Swenson (NASA ARC) and J. Mitchell (SUNY Stony Brook) for 
helpful suggestions. 


References 

1. Carr, G. C.; Erzberger, H.; and Neuman, F.: Delay Exchanges in Ar- 
rival Sequencing and Scheduling. Journal of Aircraft, vol. 36, 1999, 
pp. 785-791. 


24 



2. Bayen, A.; Tomlin, C.; Ye, Y.; and Zhang, J.: An Approximation 
Algorithm for Scheduling Aircraft with Holding Time. IEEE Con- 
ference on Decision and Control (CDC), 2003. 

3. Roy, K.; Bayen, A.; and Tomlin, C.: Polynomial Time Algorithms 
for Scheduling of Arrival Aircraft. AIAA Guidance, Navigation, and 
Control Conference and Exhibit , 2005. 

4. Bianco, L.; Dell’Olmo, P.; and Giordani, S.: Scheduling models for 
air traffic control in terminal areas. J. Scheduling , vol. 9, no. 3, 2006, 
pp. 223-253. 

5. Balakrishnan, H.; and Chandran, B.: Scheduling Aircraft Landings 
under Constrained Position Shifting. AIAA Guidance, Navigation, 
and Control Conference and Exhibit, Keystone, CO, 2006. 

6. Lee, H.; and Balakrishnan, H.: A Study of Tradeoffs in Scheduling 
Terminal- Area Operations. Proceedings of the IEEE, vol. 96, no. 12, 
dec. 2008, pp. 2081 -2095. 

7. Chipalkatty, R.; Twu, P.; Rahmani, A.; and Egerstedt, M.: Dis- 
tributed scheduling for air traffic throughput maximization during 
the terminal phase of flight. Decision and Control (CDC), 2010 49th 
IEEE Conference on, dec. 2010, pp. 1195 -1200. 

8. U.S.A.: Order JO 7110. 65T, Air Traffic Control. U.S. Dept, 

of Transportation, Federal Aviation Administration, Washington, 
D.C., 2010. 

9. Bradley, S.; Hax, A.; and Magnanti, T.: Applied Mathematical Pro- 
gramming. Addison Wesley, 1977. 

10. Vazirani, V. V.: Approximation Algorithms. Springer, Mar. 2004. 

11. Bayen, A. M.; Zhang, J.; Tomlin, C. J.; and Ye, Y.: MILP for- 
mulation and polynomial time algorithm for an aircraft scheduling 
problem. 2003. 

12. Sadovsky, A.; Davis, D.; and Isaacson, D.: Optimal routing and 
control of multiple agents moving in a transportation network and 
subject to an arrival schedule and separation constraints. In review, 
2012. 

13. Sussmann, H. J.: A maximum principle for hybrid optimal control 
problems. The 38th IEEE Conference on Decision and Control, 1999, 
pp. 425-430. 

14. Sussmann, H. J.: A nonsmooth hybrid maximum principle. The First 
European Community Nonlinear Control Network Workshop on ‘Sta- 
bility and Stabilization ’, 1999, pp. 325-354. 


25 



15. Tomlin, C.; Mitchell, I.; Bayen, A. M.; and Oishi, M.: Computa- 
tional techniques for the verification of hybrid systems. Proceedings 
of the IEEE , vol. 91, no. 7, 2003, pp. 986-1001. 

16. Branicky, M.: Introduction to hybrid systems. Birkhauser, 2005. 

17. Simic, S.; Johansson, H.; Sastry, S.; and Lygeros, J.: Towards a geo- 
metric theory of hybrid systems. Dynamics of Continuous, Discrete, 
and Impulsive Systems, Series A , 2005. 

18. Dmitruk, A.; and Kaganovich, A.: The Hybrid Maximum Principle 
is a consequence of Pontryagin Maximum Principle. Systems and 
Control Letters , vol. 57, no. 11, 2008, pp. 964 - 970. 

19. Passenberg, B.; Caines, P. E.; Sobotka, M.; Stursberg, O.; and Buss, 
M.: The minimum principle for hybrid systems with partitioned 
state space and unspecified discrete state sequence, in Proc. of the 
49th IEEE Conf. on Decision and Control , 2010. 

20. Sethian, J.; and Vladimirsky, A.: Ordered Upwinds Methods for Hy- 
brid Control. Proceedings Fifth International Conference on Hybrid 
Systems and Control, LCNS 2289 , 2002. 

21. D. R. Isaacson and J. E. Robinson III and H. Swenson and D. Den- 
ery: A Concept for Robust, High Density Terminal Air Traffic Oper- 
ations. 10th AIA A Aviation Technology, Integration, and Operations 
(ATIO) Conference, Fort Worth, TX , 2010. 

22. Swenson, H. N.; Hoang, T.; Engelland, S.; Vincent, D.; Sanders, T.; 
Sanford, B.; and Heere, K.: Design and Operational Evaluation of 
the Traffic Management Advisor at the Fort Worth. 1st USA/Europe 
Air Traffic Management R&D Seminar, Saclay, France , 1997, pp. 
17- 19. 

23. Love, D.; Horn, W.; and Fellman, L.: An evaluation of the use of the 
URET prototype tool at Indianapolis and Memphis air route traf- 
fic control centers, using benefit metrics. Digital Avionics Systems 
Conference, 2002. Proceedings. The 21st , vol. 1, 2002, pp. 1A1-1 - 
1A1-12 vol. 1 . 

24. Neuman, F.; and Erzberger, H.: Analysis of Delay Reducing and 
Fuel Saving: Sequencing and Spacing Algorithms for Arrival Traffic. 
NASA TM-103880, 1991. 

25. Sadovsky, A.; Swenson, H.; Haskell, W.; and Rakas, J.: Optimal 
Time Advance in Terminal Area Arrivals: Throughput vs. Fuel Sav- 
ings. 30th Digital Avionics Systems Conference (DASC), Seattle, 
WA, 2011. 


26 



26. Isaacson, D.; Sadovsky, A.; and Davis, D.: Scheduling for Future 
Terminal Air Traffic Operations: Problem Statement and Review of 
Prior Research. To be submitted , 2012. 

27. U.S.A.: Concept of Operations for the Next Generation Air Trans- 
portation System, Version 2.0. 2007. 

28. Papadimitriou, C. H.; and Steiglitz, K.: Combinatorial Optimiza- 
tion; Algorithms and Complexity . Dover Publications, 1998. 

29. Korn, G.; and Korn, T.: Mathematical Handbook for Scientists and 
Engineers: Definitions, Theorems, and Formulas for Reference and 
Review. McGraw-Hill, New York, 1961. 

30. Ghrist, R.; O’Kane, J. M.; and LaValle, S. M.: Computing Pareto 
Optimal Coordinations on Roadmaps. The International Journal 
of Robotics Research , vol. 24, no. 11, 2005, pp. 997-1010. URL 
http : //ijr . sagepub . com/content/24/1 1/997 . abstract. 

31. Adams, C.; and Franzosa, R.: Introduction to topology: pure and 
applied. Pearson Prentice Hall, 2008. 

32. Arnol’d, V.: Mathematical Methods of Classical Mechanics. 

Springer, 2nd ed., 1997. 

33. Jung, J. B.; and Ghrist, R.: Pareto optimal multi-robot coordination 
with acceleration constraints. Robotics and Automation, 2008. ICRA 
2008. IEEE International Conference on, may 2008, pp. 1942 -1947. 

34. Cormen, T.; Leiserson, C.; and Rivest, R.: Introduction to Algo- 
rithms. The MIT Press, Cambridge, MA, 1990. 

35. Bryson, A. E.; and Ho, Y. C.: Applied optimal control. New York: 
Blaisdell, 1969. 

36. ARINC, T. C. G.: 424-20 Navigation System Database. 2011. URL 
https : //www. arinc . com/cf /store/catalog_detail . cfm?item_id=1714. 


27 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMBNo. 0704-0188 


I 

The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection 
of information, including suggestions for reducing this burden, to Department of Defense, Washington Headquarters Services, Directorate for Information Operations and Reports 
(0704-0188), 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302. Respondents should be aware that notwithstanding any other provision of law, no person shall be 
subject to any penalty for failing to comply with a collection of information if it does not display a currently valid OMB control number. 

PLEASE DO NOT RETURN YOUR FORM TO THE ABOVE ADDRESS. 


1 . REPORT DATE (DD-MM-YYYY) 2. REPORT TYPE 

01-07-2012 Technical Memorandum 

3. DATES COVERED (From - To) 

4. TITLE AND SUBTITLE 

Efficient Computation of Separation- Compliant Speed Advisories for Air 
Traffic Arriving in Terminal Airspace 

5a. CONTRACT NUMBER 

5b. GRANT NUMBER 

5c. PROGRAM ELEMENT NUMBER 

6. AUTHOR(S) 

A. V. Sadovsky and D. Davis and D. R. Isaacson 

5d. PROJECT NUMBER 

5e. TASK NUMBER 

5f. WORK UNIT NUMBER 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

NASA Ames Research Center 
Moffett Field, California 94035-0001 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

L-jARC L-number^ 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 

10. SPONSOR/MONITOR S ACRONYM(S) 

NASA 

11. SPONSOR/MONITOR’S REPORT 
NUMBER(S) 

NASA/TM-2012-216033 


12. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified-Unlimited 
Subject Category 03 

Availability: NASA CASI (443) 757-5802 

13. SUPPLEMENTARY NOTES 

An electronic version can be found at http://ntrs.nasa.gov. 


14. ABSTRACT 

A class of problems in air traffic management asks for a scheduling algorithm that supplies the air traffic services authority 
not only with a schedule of arrivals and departures, but also with speed advisories. Since advisories must be finite, a 
scheduling algorithm must ultimately produce a finite data set, hence must either start with a purely discrete model or 
involve a discretization of a continuous one. The former choice, often preferred for intuitive clarity, naturally leads to 
mixed-integer programs, hindering proofs of correctness and computational cost bounds (crucial for real-time operations). In 
this paper, a hybrid control system is used to model air traffic scheduling, capturing both the discrete and continuous aspects. 
This framework is applied to a class of problems, called the Fully Routed Nominal Problem. We prove a number of geometric 
results on feasible schedules and use these results to formulate an algorithm that attempts to compute a collective speed 
advisory, effectively finite, and has computational cost polynomial in the number of aircraft. This work is a first step toward 

optimization and models refined with more realistic detail 

15! SUBJECT TERMS 

air transportation, arrival scheduling, separation assurance, continuous dynamics, control strategy, feasible trajectory 


16. SECURITY CLASSIFICATION OF: 

17. LIMITATION OF 
ABSTRACT 

18. NUMBER 
OF 

PAGES 

19a. NAME OF RESPONSIBLE PERSON 

a. REPORT 

b. ABSTRACT 

c. THIS PAGE 

STI Help Desk (email: help@sti.nasa.gov) 

U 

U 

u 

UU 

34 

19b. TELEPHONE NUMBER (Include area code) 

(443) 757-5802 


Standard Form 298 (Rev. 8/98) 

Prescribed by ANSI Std. Z39.18 




































